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ION  HEATING  PERPENDICULAR  TO  THE  MAGNETIC  FIELD 


ABSTRACT 

Several  theories  of  ion  heatii^  perpendicular  to  the  geomagnetic  field  are 
Imefiy  reviewed  and  assessed.  Perpendicular  heating  of  ions  leading  to  the 
formation  of  ion  conics  is  cmnmon  in  the  ionosphere  and  magnetosphere.  Ion 
conics  at  altitudes  above  a  few  thousand  kilometers  are  oftoi  associated  with 
waves  around  the  ion  gyrofiequency.  It  is  concluded  that  the  nullity  of  these 
km  ctmics  that  are  locally  heated  or  generated  over  extended  al^de  r^imes, 
may  be  best  e9q>lained  by  ion  cyclotron  resonance  heating.  At  lower  altitudes, 
particularly  in  line  r^km  of  dis^ete  auroras,  energization  by  turbulence  around 
the  lower  hybrid  fieqwmcy  seons  to  be  an  important  heating  nmchanism. 


1.  INTRODUCTION 

It  is  wdl  known  that  ions  in  the  ionoqrfiere  and  noagneto^ho-e  can  be 
heated  perpendicularly  to  the  geomagnetic  fiidd.  These  kms  may  then  move 
adiabatkaJly  up  the  fidd  lines  of  the  inhomogeneous  terrestrial  magnetic  field 
and  fixm  so>called  conics  in  vdodty  sptnoe.  These  distrfixitions  are  observed  by 
rodtets  at  attitudes  of  a  fisw  hundi^  Idltmwters  [1-7]  and  by  satdlites  up  to 
several  Earth  radii  [8-18].  The  observed  ion  ctmics  have  energies  frtmi  a  few  eV 
up  to  sevml  keV.  Mudi  of  die  heating  seems  to  occur  on,  or  near,  auroral  fidd 
tines.  It  is  now  widdy  bdieved  that  heated  ions  of  ionoqrfieric  tmgin  can 
contribute  s^nificantly  to  the  magnetoqrimic  plasma  [19]. 

Accderatitm  paraOd  to  the  geomagnetic  fidd  can  of  course  contribute  to 
the  outflow  of  ions.  The  neatly  static  imraDd  electtic  fields  that  are  bdieved  to 
accderate  auroral  electrons  downward  can  also  accderate  ions  out  into  the 
outer  magnetoqdiere.  However,  itm  distributions  with  a  distinct  conical  shape  in 
km  vdodty  space  are  t^en  olMerved,  indicating  that  perpoidicular  or  oblique 
is  an  inqxvtant  (rfienomenon.  Difl^ent  ion  heat^  medumisms  have 
been  ff»CB******  during  the  past  several  years  [20-24;  and  references  therein]. 
The  purpose  oi  flus  review  is  to  provide  a  discusdon  of  stmie  of  the  medianisms 
flurt  have  been  conaidoed  to  be  inqxntant  for  perpendicular  ion  heatmg  in  the 
terrestrial  kmoqitoe  and  magnetoqAere. 

It  a  wdl  estaUished  that  almost  all  of  the  oiergy  driving  plasma  processes 
m  the  terrestrial  nu^netosphere  comes  fi-cmi  the  sdar  wind.  This  may  be 
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conqMred  with  the  magnetosphere  of  Jupiter  where  a  large  fraction  of  the 
corre^KNiding  energy  comes  from  the  rotation  of  the  planet  [2S;  and  rrferences 
therein],  or  the  magnetos|dbere  of  a  star  whose  energy  comes  frcmi  fusitm 
processes  within  itself  Ctmsidenng  the  magnetoq^e  of  the  Earth,  energy 
originating  in  the  sohu  wind  may  eventu^  cause  oscillating  electric  fidds  at 
different  frequendes.  A  large  part  of  die  ditcussicMis  in  the  literature  ctmconing 
km  ctmics  is  really  about  what  frequency  range  of  electric  field  osdllations  is 
important  for  perpendicular  ion  heming.  Possibilities  rai^e  firom  neaiiy  static 
dectiic  fidds,  to  waves  with  fiequendes  of  the  order  of  the  ion  gyrofrequency, 
up  to  warn  above  the  lower  hybrid  frequency. 

In  the  fi:^wing  we  discuss  some  possiUe  ion  heating  medianisms.  These 
are  then  conqured  with  observatkms  in  an  attempt  to  determine  which 
mechanisms  are  more  rdevant.  One  ccmchision  is  that  cyclotron  resonance 
heating  is  inqmrtant  at  altitudes.  Furthermore,  energization  by  turbulence 
around  the  lower  hybrid  fiequency  seem  to  be  impcntant  at  low  altitudes, 
particularly  in  the  r^on  of  dscrete  auroras. 

2.  STATIC  ELECTRIC  FIELDS 

Some  dectric  fidds  may  fix  all  inactical  purposes  be  assumed  to  be 
ccmstant.  Tins  is  the  case  \dien  the  kms  traverse  ^  region  of  interest  in  much 
less  than  one  fiddcyde,  [26].  An  ion  in  a  static  homogeneous  magnetic  fidd  JBo 
wMi  an  addkkmal  perpoidicular  static  dectric  fidd  will  ExB  drift  with  an 
average  vdodty  as  discussed  m  many  textbooks.  Drifts  are 

GCHnoKm  in  the  magnetoqrfiere,  e.g.,  in  the  cusp/deR  and  in  the  polar  cap,  [27]. 
However,  sudi  drifting  distributioas  usuafly  can  not  be  confused  with  the 
essenthdly  gyrotropic  and  <^len  mote  energetic  km  ccmics. 

Locali^  dectric  fidds  with  stitmg  perpendicular  ctmoponents  are 
ormmxm  in  the  auroral  zone  [26,  2S].  Some  of  th^  fidds  might  lead  to  other 
types  of  enetgizarion  processes  than  £xB  drifts.  In  an  inhomogeneous 
perpendicular  dectric  fidd  with  diaip  enou^  gradient,  the  ion  gyromotion  mi^ 
be  disnqrted  and  the  kms  are  then  free  to  be  acederated  aax>ss  the  magnetic  field 
[22, 29].  The  nupdimim  ion  energy  adnevable  by  this  i»ocess  correqmnds  to  a 
vdodty  with  magnitude  vdtiere  .E' is  a  characteristic  value  of  the  dectric 
field.  Note  that  the  i<m  now  keeps  tins  eoargy  if  the  dectric  fidd  slowly 
disqipears,  in  contrast  to  the  £xB  drift  in  a  hcmiogeneous  dectric  fidd.  Another 
posdUlity  is  an  oblique  dectric  fidd  that  is  roughly  constam  in  a  limited  r^on 
inspoee.  The  paraOdmotkmofthe  kms  can  now  tranqxnrt  them  into  this  r^km 
dfknge  dectric  fidd.  In  crossing  the  structure  the  kms  gain  energy  equal  to  the 
potential  drop.  Ifowever,  the  increase  in  perpendicular  vdodty  is  again  limited 
by  rou^bty  tte  £xB  vdodty  obtained  from  tte  perpendicular  corrqxmmit  of  the 
dectric  fidd.  Ions  qrcruSng  less  than  one  ^rqperiod  m  the  dectric  fidd  may 
exit  with  {xtdi  uigles  near  90°  (dqrendit^  cm  the  angle  of  the  electric  fidd  with 
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re^)ttct  to  the  magnetic  Md),  while  ions  q)aKiing  iongo-  times  in  the  electric 
fidd  are  likdy  to  end  up  wi^  mcnre  fidd  aligned  vdodties  [22,  30].  Detailed 
comparisons  between  d^  and  theory  for  these  types  of  processes  are  usually 
not  possiUe  because  of  the  narrowness  of  the  regions  of  intense  dectric  fidds, 
and  the  lack  of  icm  data  with  adequate  temporal  resolution  and  sufBdent  count 
rates.  Note  that  a  static  Ej_  gives  the  same  to  all  ions,  and  thus  gives  more 
energy  to  heavier  ions. 

3.  TIME-VARYING  ELECTRIC  FIELDS  AND  CYCLOTRON 
RESONANCE 

The  sohu  wind  energy  may  directly  or  indirectly  generate  dectric  fields 
over  a  wide  range  of  fiequeodes.  Frequendes  range  fi'om  bdow  a  mHz 
(mkrqpulsations)  up  to  at  least  several  MHz  (the  plasma  fi-equency  in  the  lower 
ionoq^ere).  Many  of  these  fidds  may  be  described  as  linear  waves  in  a 
hixnogeneous  plasma  [31,  32],  but  nonlinear  effects,  sharp  gradients  and 
transient  {rfienomeoa  are  important  in  numo'ous  situations.  Possible  generation 
mechanisms  for  these  dectric  fidds  indude  both  miCToscopic  wave-partide 
interactions  (e.g.,  waves  generated  by  a  beam,  or  by  a  partide  distribution  with  a 
loss  cone  or  a  tenq)erature  anisotropy)  and  noacroscopic  phenomena  (e.g., 
Kdvin-Hdmbokz  instaUlhies).  However,  it  is  not  the  purpose  of  tins  review  to 
I»ovide  detailed  desoiptions  of  these  mechanisms. 

If  we  knew  the  electric  fidds  for  eadi  individual  ion  dong  its  orbit,  we 
could  detomine  the  "exact"  ion  distribution  by  int^rating  the  equations  of 
modem  far  each  ion.  In  practice  we  know  onfy  one  or  a  few  field  cemopements 
aloi%  the  (Mbit  of  one,  or  possfely  a  few,  q>acecaraft.  Above  the  aurord  zone,  the 
velodty  of  a  polar  orttting  satellite  is  often  roughly  perpendicular  to  Bq  and  thus 
to  the  ion  ort^.  It  is  often  assumed  that  tlm  fidds  observed  by  a  ^acecraft 
during  a  finite  time  intervd  is  represaatative  for,  but  not  the  same  as,  the  fidds 
dong  the  individud  ion  orbits.  Afthoi^  this  is  a  very  common  assumption,  it  is 
usually  rx)t  stated  e9q)licitly.  The  time  series  firom  a  qmcecraft  is  oftm  used  to 
make  a  power  qrectrum,  using  s<mie  standard  Fast  Fourier  Transform  technique 
[33, 34].  In  the  process,  details  (e.g.,  {biases)  of  the  origind  observed  signd  are 
lost,  but  this  usually  is  not  a  proUem  since  we  do  not  want  the  details  of  that 
particular  time  series  anyway.  A  serious  problon,  however,  is  that  standard 
techniques  give  informatitm  about  the  distribution  of  wave  enogy  in  fi'equoicy, 
but  not  in  wavevectcM’  space  [34].  Sometimes  wavevectors  can  be  estimate 
[35],  e.g,  by  using  a  suitable  set  of  measured  wave  firequency  spectra  together 
with  (fispersion  rdatkms  derived  from  the  observed  partide  distributions  [36, 
37].  Another  serious  problan  is  that  the  obtained  fiequoicy  spectra  nu^  contain 
s^nificant  cemtributions  firrma  essentially  static  structures  that  are  Dopplo:  shifted 
to  h^faor  frequendes  due  to  the  sa^lite  vdodty  [38,  39].  Keqmig  these 
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olmrvBtkmal  ivoUems  in  mind,  we  now  jn-oceed  to  discuss  stmie  time  vaiyii^ 
electric  fidds  that  are  inqmitant  for  ion  energizaticm. 

3.1  Resmant  Energizatitm 

A  left  drcuiarly  polarized  wave  with  fivquency  equal  to  the  ion 
gyroftequency  can  efficiently  accdente  a  positive  ion  perpendicular  to  a 
hfxnogaieous  magnetic  fidd.  Consider  an  ion  in  a  homogeneous  (infinite 
wavdength)  monochrcmiatic  left-handed  dectric  fidd  with  fi-equency/ equal  to 
the  km  gyroftequency  ff..  This  dectric  fidd  rotates  in  the  same  direction  and 
with  the  same  ai^tilar  vdochy  as  the  ion.  The  dectric  fidd  will  aj^dy  a  constant 
finrce  on  the  ion  d<mg  its  orto,  and  thus  give  a  vdochy  increase  pr(^)ortional  to 
the  time  and  hence  an  energy  increase  proportitmal  to  fi-.  This  case  is  similar 
to  die  Modenhkm  of  kms  by  a  static  dectric  fidd  paralld  to  In  this  latter 
case,  the  resonance  fiequency  can  be  said  to  be  at  zero  frequency.  Returning  to 
perpendicular  energization,  we  note  that  a  very  narrow-buMled  coherent  wave 
with  fimpiency  around  is  hard  to  obtain  even  in  a  laboratory,  and  usually  do 
not  occur  in  tte  magneto^here.  Rather,  random  phased  waves  coveriitg  a  fiurly 
tMXMd  fipequency  band  i^uding  /g  are  often  observed.  The  left-handed 
conqxment  of  these  waves  stiD  interacts  efficiently  with  the  ion.  Howeva-,  the 
km  motion  in  vdodty  r^me  can  now  be  regard  as  a  random  walk.  The 
vdodty  increase  is  prc^mitional  to  V/,  and  the  oiargy  increase  is  thus 
{nx^XHtional  to /. 

Broadband  waves  at  frequencies  around  and  bdow  the  i<m  gyrofrequency 
are  often  observed  in  the  mi^netordiere  [e.g.,  13-16,  26, 40-44].  These  waves 
are  often  associated  with  transverse  icm  heating.  Thus,  both  theory  and 
observations  seon  to  indicate  that  broadband  waves  around  fc  are  important  fin* 
ion  heatii^.  Below  we  discuss  stmae  studies  showing  that  such  cyclotron 
resonance  heating  by  waves  near  /g  is  indeed  very  impmtant  for  tnmsverse  ion 
energization.  Thb  will  be  ftdlowed  by  discustions  involvii^  possible 
contrfoutions  to  km  heating  fixxn  other  parts  of  these  broMfoand  qiectra. 

3.2.  Cyciotixm  Resonance  Heating 

One  sinq^  and  effident  itm  heating  mechanism  involves  broadband  waves 
around  the  ion  gyrofiequency  [41].  Uting  reasonable  assumptions,  Rettero’  et 
aL  [43]  apf^ed  this  cyclotron  resonaime  heating  mechanism  to  observations  in 
the  central  idasma  In  thb  shuly,  an  observed  wave  plectrum  (Figure  1) 
together  wiffi  a  Monte  Carlo  emulation  were  used  to  produce  an  ion  distribution 
quaittitativdy  nimlar  to  that  observed  by  Winnin^iam  and  Burch  [45]  at  a 
geocentric  (fistanoe  of  about  2  enth  radii  (%)  Figure  2.  (The  theory  p^cts 
that  the  kms  were  oxygen  ions.  Thb  prediction  was  confirmed  by  measuremoits 
fimn  the  mass  Electrometer  EICS  on  DE-1  [Peterson,  Khmqiar  and  Shdley, 
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private  MMumimcatioiis.])  The  thecvy  of  ion  cyclotrcm  res<»ianoe  heating  by 
teoadband  waves  is  based  on  a  diffusion  operator.  This  (^)erator  gives  tbs 
dffiision  rate  correspcmding  to  the  rand(Mn  walk  of  ions  in  vdo^  space.  Usii% 
the  kmg  wavelength  api»x>ximati<»i,  the  op«ator  can  be  written  in  a  very  ample 
fimn.  This  ^^yroximarion  shcHild  be  val^  for  examfde,  for  tvoatfoand  Alfvdi 
waves  in  most  r^kms  of  the  magnetosphere.  The  approimate  diffiision 
coefficient  is  then  sin^yly  proptvtional  to  the  dectric  fidd  spectral  density  (S)  at 
the  local  gyrofiequency  of  tte  ion  specie  of  interest  [41,  43].  The  resulting 
average  heating  rate  per  ion  (Q)  can  also  be  shown  to  be  fn^portional  to  the 
iq)ectral  density  (and  rluM  to  the  noean  square  of  the  dectric  ^d  fluctuations)  at 
tte  local  ion  gyrofiequency  [41]; 

Q-(q2/2m)SL  (D 


Hgure  1.  TyfHcal  electric  fidd  spectral  doiaty  in  the  coitral  plaana 
sheet.  Geocentric  distance  »  2  R^,  invariant  latitude  »  60*’.  (M. 
Mdlott  and  D.  Gumett,  private  communication,  1986) 
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9  rad  m  affe  the  charge  and  mass  of  the  imi,  respectively.  The  qiectral 
drasity,  Sl,  in  etpiatkm  (1)  is  the  fraction  of  the  spectral  density  diw  to  h^-hand 
polari^  waves.  As  an  example  of  the  rdative  effictoicy  of  this  heating 
mechanism,  we  note  friat  a  q>ectral  density  of  10  (mV/m)2/Hz  would  yidd  an 
average  heating  rate  of  30  eV/s  for  oxygen  ions  [13].  This  heating  rate  is  high 
enough  to  explain  km  distributions  at  essratially  90  degrees  pitch  angle  in  the 
dayskte  magnetosphere  at  an  altitude  of  about  13,500  km  (Figure  3  ai^  4). 
Figure  Sa  dqpicts  how  an  observed  wave  spectrum  may  be  used  to  estimate  the 
local  ion  hea^  rate. 


Figure  2.  Bottom  panel  is  a  contour  diagram  of  the  observed  ion 
conic  distiibutitm  flmction,  measured  by  the  HAPI  instrument  on  the 
D£-l  satellite  (Wimrin^fiam  and  Bur<^  1984,  [45];  the  fiun  that  the 
ions  were  oxygra  ions  were  confirmed  by  the  EICS  instrument  on 
DE-1,  Peterstm,  Klun|)ar  and  Shdley,  private  communications). 
Top  pand  is  the  calculi^  ion-vdocHy  dikribution  (Rettoer  et  al., 
1987,  [43])  uang  the  dectromagnetic  ion  cyclotron  resonance  theory 
(Cha^  et  al.,  1986,  [41])  plotted  in  the  same  way  as  the  observed 
distribution. 
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V^ao^m/s)  VjOO^nj/s) 

Figure  3.  Distribution  functi(His  of  heated  ions  observed  by  the  Viking  satdlite 
in  the  cu^cleft  re^^n  (Andr6  et  al.,  1988,  [13]).  The  left  distribution  is 
obtained  near  the  equattn-ward  edge  of  the  cusp/cleft  (2231:14  to  2231:34  UT 
(Ml  July  24,  1986)  ajnd  is  measured  simultaneously  with  the  intense  broadband 
waves  in  Figure  4. 


Figure  4.  Wave  dectric  fidd  spectra  obsmved  by  the  Viking  spacecraft  in  the 
cusp/deft  r^km  (Andr6  et  al.,  1988,  [13]).  The  more  imense  waves  (2231:27 
UT)  are  observed  simultaneousty  with  the  heated  ions  in  Figure  3. 


Q0CS1S2 
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Q«Sn^ 


RgureS.  Skttdies  delating  how  an  observed  wave  spcctnim  can  be 
used  to  esthnate  local  ion  heatiug  by  the  (a)  ^dotron  resonance,  (b) 
douMe  c^dotron  resonance,  and  (c)  sutdiannonic  resonance 
medianisais.  Shaded  areas  depict  bandwidths. 
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Actually  the  rescmance  condition  lidf  =  AijiVji  should  be  used  to  find 
out  at  which  fi-equendes  if)  ions  are  in  resonance  wHn  a  wave.  Here  A:||  and  V|| 
are  the  parallel  wavevector  and  the  parallel  componoit  of  the  particle  vdodty, 
vdule  JV  is  an  integer.  This  relation  takes  into  account,  e.g.,  the  finite 
wavdo^lth,  and  the  posabilhy  of  interactions  at  higher  harmonics.  For 
emissiQns  with  long  wavctengths  sudi  as  most  waves,  the  effect  of  the 
term  i:||V||  becomes  small.  Interaction  at  higher  harmonics  often  require  a 
nonzero  perpendicular  component  of  the  wavevector. 

3.3  An  Exanq)le:  Cyclotron  Resonance  Heating  in  the  Cusp 

As  the  ions  move  up  the  fidd  line,  they  interact  with  the  spectral  dendties 
of  different  fiequency  bai^  as  the  local  ion  gyrofi'equency  dianges  due  to  the 
inhomc^eneity  of  tte  geomagnetic  fidd.  As  an  example  of  a  test  of  the 
cyclotron  resonance  heating  medbanism,  we  discuss  some  observations  in  the 
cusp/deft  r^on  of  tl»  magneto^here.  These  studies  include  nuq>ping  of 
q)ectral  densities  down  the  geomagnetic  fidd,  and  estimates  of  the  ion  hotting  at 
various  altitudes  as  these  partides  move  upward  [13, 14, 46].  The  cusp/cleft  is  a 
r^(m  of  the  dayade  high-latitude  magnetoq)here,  extending  fiom  less  than 
0900  to  more  than  1500  magnetic  local  time  and  with  a  latitudinal  width  of  a  fisw 
degrees  [47-49].  In  this  r^on,  predi»tation  of  shocked  solar  wind 
(magnetosheath)  plasma  usually  occurs.  Low-firequency  broadband  waves  are 
also  crmumm  in  this  r^;ion  [SO,  51].  A  sketch  showing  a  satdlite  moving 
poleward  into  the  cuiq>/cleft  region  is  shown  in  figure  6.  Here  region  A  is  the 


Figure  6.  Sketdi  showii^  upgoii^  iras  (solid  arrows)  moving  into 
the  cuq>/deft  (r^on  A)  fix>m  the  equatorward  region  (B)  due  to 
convectkm.  Shading  in  r^km  A  indicates  low-frequency  waves 
whidi  can  heat  the  ions.  Fr<»n  Andr6  et  d.,  1990,  [14]. 
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cu^deft  with  rather  intense  broadband  low-frequency  waves  (shaded),  and  also 
with  injected  magnetosheath  particles  (got  indicated).  In  the  equatmnvard  r^on 
B  the  wave  intensity  is  much  lower.  Upgmng  ions  are  shown  by  solid  arrows. 
Since  there  ofren  is  a  high-latitude  poleward  convection  field  [52],  these 
ions  can  drift  firom  B  into  r^ion  A.  A  poleward  moving  satdUhe  may 
observe  upflowing  low-energy  irms  at  point  1  in  Figure  6.  Waves  and 
crmvectkm  drift  can  be  detected  at  point  2  and  may  be  m^>ped  down  to  estimate 
ion  heating  and  drift  at  lower  altitudes.  The  heated  ions  can  be  obsoved  at  point 
3. 


I%ire  7.  Data  fi:om  an  evat  where  DE  1  moved  poleward  into  the 
cu^cleft.  The  top  pand  shows  dectric  fidd  q)ectral  density  in  units 
of  V^m’^Hrl  (the  O*  gyrofiiequenr^  is  0.85  Hz).  The  second  and 
third  pands  show  count  rates  and  O*  ions  with  energies  (E) 
from  10  eV  to  17  keV.  The  left  and  right  scales  on  the  gr^  bar 
correqxmd  to  die  observed  H**  and  CX*'  count  rate,  req)ectivdy,  in 
counts  pCT  sanq^,  rdsdi  is  proportional  to  number  flux.  Thelfourth 
pand  ii^cates  the  (Htdi  angle.  The  bottom  pand  shows  count  rate 
(again  in  counts  per  sanq)le)  of  imis  with  energies  bdow  65  eV  as 
a  fimction  of  qrin  angle  and  time,  'indiere  the  solid  line  indicates  the 
upgoing  direction.  Note  the  onset  of  waves  (first  pand)  and 
energetic  (X**  kms  (third  pand)  at  about  0124  UT,  and  the 
disiypeatance  of  low-eietgy  (X^  kms  (lowest  pand)  at  the  same  time. 
Frmn  Andrd  et  d.,  1990,  [14]. 
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Rgure  7  shows  data  from  DE>1  obtained  when  the  satellite  is  moving 
pdeward  into  the  cu^cleft,  i.e.,  simUar  to  Figure  6.  The  top  pand  in  Figure  7 
shows  the  low-frequency  dectric  field  q)ectral  density.  The  wave  mtoisity 
mcreased  after  abom  0124  UT.  In  particular,  the  qwctral  density  at  the  O'*' 
gyrofinequency  (0.85  Hz)  increased,  and  this  is  imp(»tant  for  the  resmiant  heating 
we  cmidder.  At  about  0124  UT  the  satdlite  entered  the  cu^/deft  r^ion,  and 
this  corresponds  to  the  edge  between  r^ion  A  and  region  B  in  Figure  6.  The 
second  and  third  panels  of  Figure  7  show  magy  time  spectrograms  for  hyd  ^ 
and  ox^en.  Insi^  the  cusp/deft,  the  characteristic  hydrogen  "butterfly  pat  a 
caused  by  the  poleward  ExB  drift  of  ii^ected  ions  can  be  seen  [52].  Intense, 
rdativdy  local  heatii^  of  upflowing  ions  is  indicated  by  the  cmiic-type 
distributions  with  a  density  of  roughly  1  cm~3  that  q>pear  in  ibc  third  pand  at 
about  0124  UT.  It  is  the  heating  of  these  ions  we  consider  in  some  det^.  The 
bottom  pand  of  Rgure  7  shows  the  counting  rate  of  ions  with  energies  above 
the  ^Moecraft  potential  and  bdow  about  65  eV.  The  data  are  presented  in  an 
angie*time  ^MCtrogram  format,  and  iixiicate  that  upflowing  oxygen  ions  are 
present  before  about  0124  UT.  After  this  time,  no  low-enogy  O*"  irnis  are 
observed.  Cmrqraiison  with  the  energetic  m^en  spectra  in  the  third  pand 
suggests  that  the  cool  O'**  ions  are  atergaed  out  of  the  low  energy  ngon  (bdow 
65  eV)  after  about  0124  UT,  i.e.,  inside  the  cuqi/deft. 

A  Monte  Carlo  simulation  has  been  perftwmed  to  test  the  idea  of  cydotrmi 
resonance  hearing  in  detail  [14].  The  dectric  Adds  causing  pdeward  convection 
and  ion  heating  were  observed  along  the  satdlite  tradr.  Usitig  reasonatde 
assutrqitions^  these  dectric  Adds  could  be  mapped  down  the  gemnagnetic  Add. 
Oxygen  kms  were  then  laundied  at  dififertent  altitudes  along  the  equatorward 
ec^  of  the  cuqr/deft.  Eadi  partide  was  Allowed  in  small  time  stqrs  up  to  the 
sateiltte  altitude.  Tor  each  time  stq>  the  poleward  drift  due  to  convection, 
the  upward  motion  due  to  the  patelld  \docity  and  the  adid)atic  Adding  were 
calculated.  In  additimi,  the  dumge  in  perpendicular  vdodty  due  to  interaction 
with  waves  around  the  local  gyrofrequency  were  estimated  by  iismg  the 
Monte  Carlo  technique. 

Figure  8  shows  the  wave  input  parameters  to  the  siiiiulation  together  with 
the  energies  of  the  observed  and  caloilated  O'''  distributions.  The  equatorward 
edge  ofthe  cu^deft  is  at  about  0124  UT.  The  middle  pand  shows  the  dectric 
Add  qMctral  dmsities  at  the  O^  gyrofrequency  fi>r  eadi  satdlite  qnn  together 
with  mean  values  to  guide  the  eye.  Points  equatorward  of  the  edge  are  included 
for  cmnparison.  The  lower  pand  diows  the  ^>ectrd  indices  a  (where  S  oc /*<>) 
used  to  nuq>  tiie  wave  intensity  to  Iowa-  altitudes.  The  upper  pand  shows  the 
mean  energy  correqxxxliiig  to  the  obsoved  hot  O^  distri^ons  (solid  drdes) 
and  from  the  sinadatimi  (open).  The  two  curves  are  in  good  agreement  up  to 
dxwt  0130  UT.  After  0130  UT,  the  east-west  drift  increases  significantly, 
indicating  a  turmiig  cmivection  pattoiL  Thus  the  problem  becmnes  three- 
dimensional  (particles  drifting  not  only  poleward  but  also  eastward),  vriiidi  is 
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one  reason  the  agreement  between  theory  and  oteervations  is  not  so  good 
at  lagher  latiti^. 


UT(s)  4800  SOOO  5800  5400  5600 

0T(fa3iu)l£80O0  1.83:20  1:85:40  1:30:00  1:3320 

z(km)  >517  0  517  1033  1550 


Kgure  8.  The  iqiper  pand  shows  cusp/deft  mean  omgies 
obtained  by  DE  1  <ni  ^gust  10,  1984  (sdid  circles),  and  frmn  a 
sinKdatkMi  (<^)en  drdes)  as  a  ftinction  time  (UT  given  bodi  as 
seocmds  aftCT  midniglit  and  as  hours  and  miniites).  The  middle  pand 
shows  the  dectric  fidd  ^lectral  density  5  at  die  Q*'  gyrofiequency. 
hfom  values  are  also  indeed  in  tins  pand.  The  lowest  pand  shows 
the  qiectnl  index  a.  The  satdlite  moves  pdeward  into  tte  cxup/ckA 
at  about  0124  UT  vdiere  the  wave  intensity  mcreases  and  ion  heathy 
starts.  From  Andrd  et  d.,  1990,  [14]. 

Observed  (fistributkm  fimctions  are  diqdayed  in  the  uppa  pand  of 
Rgure  9.  The  left  distribution  is  obtained  near  the  equatorward  edge  of  the 
cmp/deft,  and  seems  to  be  rather  locally  heated.  Ihe  right  distributkm  is 
observed  Inside  the  cu^deft  and  indicates  name  nrodocd  heating.  Ikretbeions 
are  more  afifocted  by  adisbaric  fdding,  and  they  are  more  "fifted”  in  energy 
toward  higher  paraUd  vdodties.  This  is  what  we  expect  since  ftiese  partides 
have  qpent  more  time  in  the  cuqj/deft  waves  and  have  traversed  a  grei^ 
aUtude  m^.  The  lower  part  of  Figure  9  shows  the  correspond  distribution 
nyww  411*  Miwilarinn  DetsSed  omtqMuison  thew  (hstrftxitions  wifti 
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tile  obsefvatkms  is  not  meaningful  since  the  energy  diannds  of  the  mass 
Electrometer  are  rather  wide.  Howevo’,  the  sinmlation  distributimis  show  the 
same  tendency  as  the  observatimis  in  the  sense  that  the  equatmward  distribution 
is  more  locally  heated.  Thus  there  is  good  agreement  between  observed  O'*' 
mean  oiergies  and  distribution  functions,  and  the  corresfXHxling  prc^ierties 
obtained  from  simulatitms  using  observed  parameters  as  an  ii^. 
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Figure  9.  Distribution  functions  of  O'*'  ims  observed  by  the  EICS 
mass  E>«ctrometer  on  DE-1  (ui^ia'  pands)  and  obtained  by  a  Monte 
Carlo  simulation  (lower  pan^).  Itiere  are  two  ctmtours  of  [diase 
qiace  density  per  deca^  and  the  innermost  (highest)  contour 
ccKTeExmds  to  2  x  10^  s^/knfi.  The  density  of  the  simulatkm 
distijbutkms  is  taken  to  be  1  cm'^,  vdudi  is  the  ai^iroximate  density 
oftiie  observed  distributions.  The  left  distributions  are  obtained  fixmi 
0124:45  to  0126:21  UT,  i.e.,  near  the  equatorward  e^  (compare 
with  Figure  8),  while  the  ri^  distributions  are  from  0128:45  to 
0130:21  UT.  Note  that  the  equatorward  distributions  are  mart 
locally  heated  than  the  distrfiiutions  from  inside  the  cuE)/deft.  From 
Andrdetal.,  1990,  [14]. 
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fo  the  above  discusskn  of  kn  heating  in  the  cusp/deft  we  assume  that  CX'' 
ions  witfi  energies  of  a  few  eV  exist  at  altitudes  of  a  few  thousand  kilonitters. 
These  ions  are  then  heated  to  of  a  few  hundred  eV  on  their  wi^  iq)  to 

10,000  or  20,000  kflometers  altitude.  Hus  assumption  is  conastent  with 
Goqugate  observations  of  ions  by  the  Akdxmo  satdlite  («  5,000  km)  and  the 
OE^l  spacecraft  («  20,000  km)  [S3].  Other  observatkms  by  single  satellites 
otMD&m  fee  occurrence  of  ion  heating  up  to  a  few  eV  in  fee  low  ahitutfe 
oap/ckA  (see,  e.g.,  ion  data  below  10  eV  from  the  DE-2  satellite  at «  600  km  in 
nate  2  of  Majnterd  et  al.  [54]). 

3.4  More  About  Cydotron  Resonance  Heating 

Carefid  tests  of  the  cydotnm  resonance  heating  medtanism,  involving 
simultaneously  observed  wave  and  k»  data,  have  been  perftnmed  also  in  the 
n^itside  aur^  zone  [44].  Ifere  drifts  can  often  n^^ed,  and  now 
fee  hemmg  piYtcess  may  often  be  diaracterized  1^  two  parameters.  These 
parameters  correspond  rougUy  to  the  velocity  a^  pitdi  an^  of  the  ion 
featiiNition,  and  may  be  independendy  detennined  from  the  wave  and  paitide 
observations  of  a  conic  event.  The  study  Crew  et  al.  [44]  again  shows  good 
agreement  between  observations  and  cydotm  resonance  heating  theoiy. 

Statistical  studies  of  kn  mass  spectrometer  data  fimn  both  DE>1  [17]  and 
Akdxmo  [18]  indicate  feat  kms  are  not  suddenly  heated  at  (me  altitude  and  then 
adiabaticall^  move  upwvd.  Rather,  the  pittfe  an^  of  the  km  chstributkms 
suggest  heightHotegiited  transverse  accderatkm  of  kms  over  a  wide  altitude 
range.  Again  this  is  consistent  wife  the  scenario  of  cydotrcm  resonance  heating 
presented  above. 

We  note  that  no  spectral  peak  around  ^  is  needed  frxr  cydotrtm  rescmance 
heating  to  be  effident.  Obviously  there  must  be  a  siginficant  amount  of  left- 
handed  waves  around  fee  gyrafrequency.  Often  an  assumed  fiactxm  of  about 
10%  of  leftfeanded  waves  is  suffic^oit  to  exjriain  the  observed  ion  heating  [14, 
41].  The  exact  wave  polarizatkm  is  hard  to  determine  from  electric  fidd 
observations  but  the  hbove  assunytkms  mem  ccmsistent  wife  magnetcnneter  data 
[14].  Furthermore,  we  note  that  usualfy  no  decrease  in  qpectral  density  around 
the  local  is  observed.  Such  a  decrmue  mi^  at  first  be  eiqpected  since  the 
waves  are  transArii^  eneigy  to  the  kms.  However,  it  must  be  realized  feat  only 
a  mmor  part  of  the  observed  wave  intensity  o(mtributes  to  km  eneigizatkm 
(exdudiag  right-handed  waves  and  electric  fi^  Doppkf  shifted  due  to  satellite 
molkm).  Fbrfeermore,  fee  observed  partide  oiergies  are  often  obtained  by 
continuous  heating  wto  fee  ions  move  thousands  of  kfiometers  iq>  the 
geomagnetic  field,  corresponfei^  to  a  wide  rai^(rf‘gyn^equeocies.  Tlius,itis 
not  stoptisteg  feat  obsenutkms  often  feow  just  a  broadband  wave  q)e(mfum 
associated  wife  km  conics. 
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Most  (rf*  the  observttioas  discussed  above  are  frcHn  altitudes  of  a  few 
dxMisand  iq>  to  about  20,000  km.  At  first  one  might  esqiect  the  heatmg  rate 
to  aim|dy  increase  with  altitude  (decreasing  gyrofi«quency  and  this  increasii^ 
spectral  density  in  the  typical  broacfiiand  qiectrum).  At  some  stage, 
approodmaticwM  used  in  the  theory  will  of  course  break  down.  For  examfrie,  the 
gyroradius  becomes  so  large  that  the  ^cm  can  not  be  assumed  to  be  locally 
homogeneous.  FurtheroKMe,  there  are  indicmions  that  the  Alfedn  wave  spectrum 
has  a  peak  at  smne  small  but  nonzero  firequency.  Maxiimim  heating  then  occurs 
where  this  peak  matches  the  local  gyrofirequency.  Indeed,  observations  with 
good  fi:e9iency  resolution  made  by  the  Viking  sateOite  in  the  cu^cleft  [55]  and 
in  die  auroral  zone  [16, 26],  and  by  the  Intercosmos-Bulgaria  1300  spacocnft  in 
the  aurcaal  zone  [56],  indicate  sudb  a  ^lectral  peak  between  0.1  and  1  Hz.  Hus 
corresponds  to  a  maximum  heating  rate  of  CX*’  kms  at  an  altitude  of  a  few/^. 

At  altitudes  bdow  a  few  thousand  kilometers,  cyclotrai  rescmance  heating 
may  be  less  inqxxrtant  sinqily  because  the  qiectral  demity  at  the  local  often  is 
rafter  low  (the  gyrafiequency  is  high)-  Thus,  we  cmiclude  that  cydotrcm 
resonance  heating  is  impntam  at  fthudes  firom  a  few  thousand  kilometors,  and 
at  least  tqi  to  a  fisw  R^.  Later  we  discuss  some  mechanisms  that  may  be 
inqiortant  at  low  ftitudes.  Hrst,  however,  we  discuss  the  posable  contrftution 
to  ion  heating  by  a  broadband  spectrum  fixnn  ofter  fiequendes  than  the 
^it^equency. 

4.0  HEATING  BY  WAVES  BELOW  THE  ION  GYROFREQUENCY 
4.1  Multqiie  Cydotrai  Resonance  ifeating 

It  has  been  atggested  by  Temerin  and  Roth  [42]  that  nmUnear  resonance 
could  also  produce  transverse  ion  heating.  Oneprooessinvdvestwo  waves  with 
finite  wavoaimbers,  vdiere  the  sun  of  die  fiequacies  matdi  the  gyroftecpiency. 
This  mechanism  may  be  called  "douWe  cydotron  abscuption”  (or  ”do(d^ 
cyclotron  resonance*)-  Assumng  low  anqplitude  dectiic  fluctuatkms,  the 
average  heating  rate  per  ion  fiir  dus  process  can  be  shown  to  be  proportion^  to 
dm  product  of  die  dectric  field  spectral  densities  at  the  two  fiequendes. 
Aldim^  dus  mechanism  is  less  effideot  than  ion  cydotron  resonance  heating 
(since  it  is  of  hi^ier  order  in  spectral  demities),  it  can  involve  the  entire 
fiequenQT  range  (bdow  the  local  gyrofiequency)  of  the  spectrum  and  can 
produce  a  fiur  amount  of  transverse  ion  acoderation.  Figure  5b. 

Assuming  that  waves  within  smne  bandwidth  around  half  the  i<m 
gynfteqpienqr  oooitfeute  to  douUe  cydotnm  absOTption,  Ball  [57,  58]  obtained 
an  expression  fijr  the  average  heating  rate  per  ion; 
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Here  /3»  2icy^  ^  is  the  rehoive  bandwidth  around  half  the  gyrofrequency 
over  wfaidi  doufaie<ycioCron  absorptkm  is  assumed  to  occur,  kj^  is  the 
perpeodkubr  wave  numbo*  of  the  waves,  is  the  dectrk  ^>ectral  density 
around  half  the  gyrofiequeacy,  «ad  rjocA  ^  *  dimensioiiless  &i^or  ^1.  The 
spectral  denshy  is  assumed  to  be  ti>e  same  aU  fiequendes  of  interest,  and  thus 
tte  spectral  density  squared  occurs  in  Eq.  (2).  The  derivation  of  rdation  (2) 
reqinres  tlMt  kx  be  mu^  less  than  the  inverse  of  the  ion  gyr<Madius  and  that  o^ 
waves  nev/^  ccmtfftute  to  the  heathy.  These  assunqitions  seem  necessary  to 
obtain  a  u»^  csqMWSskm  fin  the  double-cyclotron  heating  rate,  but  are  not 
essential  features  of  this  heating  process.  The  fiictor  is  included  in 
expressioo  (2)  because  it  is  possfele  feat  not  aU  the  spectral  density  observed  at 
f/2  cmttrfeutes  to  this  heafeig.  As  in  the  case  of  c^otitm  resonance  heating, 
for  exanqd^  Dopfder  shift  due  to  the  spacecraft  motion  of  qurtial  dectric  ^d 
structures  nu^  occur.  We  note  that  double-^clotrmi  absenption  is  very 
ineffideot  fin  long  wavdengfes  (small  i^x)-  perpendicular  wavevecten  is 
often  very  lund  to  estimate  fiom  observafems.  However,  an  iq>per  Ihnit  of  the 
heatmg  rate  given  by  (2)  can  be  obtained  by  using  a  ctnre^pmiding  to  the 
inverse  trf*  the  ion  gyroradhis. 

The  estimate  of  the  double-cydotrcm  absorption  heatiqg  rate  in  Eq.  (2)  is  a 
feettn  of  S  larger  than  fee  result  given  by  Temerin  and  Rofe  [42].  A  possiUe 
ex|danation  for  this  discrepancy  is  discussed  by  Bafl  [ST].  However,  uncertainties 
in  fee  estimates  obtained  from  observatkms  c€  parameters  such  oakj^  ate  rather 
large.  Thus  the  discrquuicy  between  Temerin  and  Roth  [42]  and  Ball  [57,  58]  is 
not  sufiBdently  huge  to  influence  an  assessment  of  the  possMe  inqxHtance  of 
doubie-cydotron  absorption  as  a  source  crmics. 

A  detailed  oonqNuison  of  douUe^ydotrcm  absmption  wife  cydotnui 
resrmance  heatirig  in  tte  oasp/ckA  r^on  of  the  magnetoqfeere  was  perfixmed 
by  Ball  and  Andrd  [59].  Hus  study  indiuled  the  DE>1  data  in  Hgures  7  and  8, 
a^  used  the  heating  rates  fuesented  by  Ball  [57,  58].  It  was  conduded  that 
douUe-cydotrcNi  drsorption  migbt  give  an  appredaUe  contribution  to 
energiatimi,  e^redally  vdien  tiie  heating  was  accomplished  loealfy  (in  altiturfe). 
However,  ^  nu^  ermdusion  of  Ball  and  Andr^  [59]  is  still  that  cydotron 
tesrmance  heating  is  the  more  effident  medumism  for  the  oxygen  conics 
ol»erved  in  the  oasp/deR.  r^otL 

Fnxn  a  thecxdical  point  of  view,  th^e  exist  obviously  other  higher  ordha 
rescmance  fxocesses  invdvir^  spectral  densities  at  nx>re  tl^  two  fiequendes. 
However,  tile  heating  rates  under  these  drcumstances  will  be  of  even  higber 
order  in  wave  anqditudes  (or  ^lectral  douities).  For  reasmiably  low  wave 
ampfitudes;  tiisse  hitler  ordv  eflbcts  will  be  much  less  effident  and  can  usually 
ben^ected. 
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4.2  SuUianiioiiic  Heatii^ 

From  die  discusskm  in  the  fvevious  sectkMi,  h  follows  that  waves  with 
fifetpiencks  /  around  the  Mh  onto  sifohamiomc  if  -  ft/N,  N  ^  2,3,4...)  can 
resonant^  heat  ions  for  fimte  wavdength  dectric  fidd  fluctuatkms.  Figure  5c 
[42]. 

The  average  rate  is  then  propottionai  to  the  Mh  power  of  the 

dectric  fidd  spectral  density,  5^,  at  tte  siddiamxMaic.  For  sufficiently  low 
amiditu^  fluctuating  fidds,  this  indicates  the  fdlowing  c<mdusi(His;  (1) 
Cydrnrcm  resonance  (M  =  1)  is,  fix  tyjncal  broatfiiand  spectra,  often  more 
efficient  than  subharmonic  resonances,  (2)  2nd  orda-  subhamKMuc  restmance 
heating  (vdudi  is  equivalent  to  the  speciai  case  of  douUe  cydotrmi  resonance 
fix’  /]  » fJ2)  is  generally  more  efi^ent  than  those  due  to  other  hi^ier  order 

subharmonic  resonances,  and  (3)  For  Ivoadband  dectric  fidd  spectra,  it  is  mmre 
inqxHtant  to  indude  all  die  douUe  cydotixm  resonance  effiects  (in  addhimi  to  the 
N  ^  2  au^bumoBic  tfkci)  than  to  include  other  hi^ier  wder  (N  >  2) 
suMiannonic  effects. 

We  end  this  section  by  noting  that  recoitiy  Terasawa  and  Nambu  [60]  have 
crasidered  suMiarmonic  km  heating  by  magnetosonic  waves.  However,  this 
process  requires  magnetic  fluctuadtms  larger  than  those  usually  detected  in  the 
aurmal  and  cuq[>  rqkms. 

4.3  Nonresonant  Enogization 

Waves  with  fiequendes  much  bdow  the  km  gyroflequeoc^  may  eneigize 
kms.  Since  bodi  die  suUutrmonic  heating  and  multi|de  cydotron  restmance 
heating  wifi  be  of  very  high  order  in  wsw  amplitude  here,  tfa^  can  generally  be 
nqglected.  It  has  ben  si^gested  that  such  waves  instead  m^  nonresonuitly 
cause  significant  km  heating  in  the  magneto^here  [16, 61]. 

To  test  this  idea  of  nonresonant  energizatkm.  Ball  and  Andrd  [59] 
peifonned  a  numoical  calculation  based  on  a  modd  of  an  observed  broadband 
wave  plectrum.  The  ufqier  pait  of  Flgine  10  dqncts  sudi  a  modd  iqiectrum, 
contnicted  by  siqierinqiosing  many  sine  waves.  A  qpectium  similar  to  this  was 
observed  by  the  Viking  satdhte  in  ctn^mctkm  with  keV  protons  [15],  The 
modd  waves  are  lineaiiy  pdarized,  and  hence  indude  a  left-handed  component. 
It  was  assumed  in  this  ^culation  that  the  wavdengths  of  the  fluctuations  w^e 
veiy  kmg,  *  0,  so  that  muhqile  cydotron  resonance  and  subharmonic 
resonances  are  not  present. 

To  test  the  idea  c(  nonrestmant  heatup  the  model  spectrum  was 
subffivided  inu>  three  parts;  (1)  <me  resonant  part  (RES)  around  the  protcm 
gyrofirequency  at  30.0  Hz,  (2)  one  low  fi-equoicy  (LF)  nonresonant  part,  and  (3) 
(me  h^  fie^iency  (HF)  mmrescnumt  part.  (Thus,  the  nonresonant  part  LF 
contdnsmoreflianjusttte  very  low  fiiequency  part.)  The  efifoct  of  die  differoit 
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parts  of  liw  q>ecmim  on  icm  eoa^izatkm  was  tested  by  stnq^  striving  the 
ecpiation  of  motion  fw  (wotons  subject  to  the  (known)  time-varying  dectiic  fidd 
of  the  modd  in  the  |m»enoe  of  a  static  uniform  magnetic  Md.  The  resuh 
obtained  by  averaging  over  many  ions  is  di;q;>]ayed  in  the  lower  panel  of  Figure 
10,  diowmg  the  mean  proton  mergy  versus  time.  The  diffidence  between 
restmant  heatii^  by  waves  around  the  proton  gyrofrequency  and  nonrestmant 
energizatkm  by  other  parts  oi  the  q)ectrum  is  dramatically  Olustrated.  It  is  seen 
that  the  low  frequency  part  of  the  qjectrum  causes  the  perpendicular  {Moton 
energy  to  increase  nqridiy  to  around  2S  eV,  after  which  no  further  net  energy 
transfer  occurs.  The  energization  {ffocess  due  to  the  frequency  part  of  the 
spe^rum  is  similar.  The  resuhiqg  energy  is  lower,  however,  since  the  q>ectral 
density  is  lower  there.  Estimates  show  that  including  wave  power  at  even  higher 
fretpiaoctes  will  not  fnovide  much  a^htkmal  km  energization. 


Rgure  10.  Proton  energizati(m  due  to  dectiic  fidd  fluctuations 
condstent  with  a  spectrum  obsoved  togedimr  with  keV  (Hoton 
conics.  The  upper  pand  shows  the  qmctral  density  cdculated  from 
the  dectric  Adds  u^  in  the  numeritdl  calculations.  The  qiectrum 
has  been  qriit  into  three  parts,  a  resonant  part  of  width  f/S,  and 
nomeaomart  {Wts  at  lower  and  higher  (iq>  to  2fc)  frequences.  The 
lower  pud  dxnvs  the  average  energy  cf  protons  subjected  to  the 
dtoctric  fidd  ffiictuati<ms  in  these  frequency  r^ons.  From  Ball  and 
AM,  1991,  [24]. 
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The  iverage  km  emrgy  <A>tained  due  to  fluctuations  in  a  frequency 
r^(ion  well  removed  from  i.e.,  outside  roi^My  0.6  fc  </<  1 .4/^,  can  easily 
be  [59].  A  roi^  cakuiatkm  inclu<k»  a  static  perpendicular  dectric 

Add  of  anqplitude  equal  to  the  intagrated  nonrestmant  density,  and  thus 

Ql<-(2lii/B^j^SUW  (3) 

Here  is  the  absdute  value  of  the  ambient  magnetic  Add.  This  situation  may 
be  ccmipared  wifli  a  linear  osdllatCH',  which  is  externally  driven  at  a  frequency  fru* 
awi^  res(mance.  This  osdlhOor  will  reach  a  maximum  amplitude,  ii^ch 
then  does  not  increase  with  time. 

Returning  to  Figure  10,  we  note  that  after  a  few  gyroperiods  the  mean 
prot<m  eneigy  due  to  resonant  heatii^  by  waves  around  fc  is  higher  than  the  total 
enogy  obtained  fiom  nonrescmant  mo-gizaticm  processes.  Again  we  note  that 
no  peak  in  the  jq)ectnun  around  the  gyrofirequency  is  needed  for  effident 
res<mant  heating,  {xrovided  suflfident  wave  power  is  present  near  the 
gyrofiequeocy. 

The  numerically  obtained  heating  rate  agrees  well  with  the  theoretical  value 
from  Eq.  (1),  as  discussed  by  Ball  and  Andrd  [59].  Sudi  restmant  heating  can 
readily  expl^  the  observed  keV  protons  reported  by  Ihiltqvist  et  al.  [15].  On 
the  oUier  Hand^  mmresonant  cneigizaticm  processes  can  cmly  account  for  a  small 
fraction  of  the  observed  km  energy.  We  ccmdude  tto  fundamental  ion 
cyclotron  resonance  heating  can  easily  exfriain  observed  hi^  energy  proton 
conics  (above  100  eV)  whidli  are  hard  to  esq^hm  with  ai^  other  mechanism. 

For  events  involving  heavior  ions,  the  atuation  is  less  dear.  As  discussed 
above,  detailed  studies  of  heating  events  show  fliat  resonant  heating  by 
obsenmd  broadband  waves  can  heat  the  ions  to  the  observed  energies. 
However,  since  the  gyrofrequendes  of  heavy  kms  are  raflmr  low,  it  will  be  more 
difiBcult  to  dififiaenttate  the  contributions  of  individual  energizatitm  mechanisms. 
At  least  in  some  instances,  double-cydotron  absorption  may  significantly 
contribute  to  CX*'  heatup  [59].  Similaily,  ncmresonant  intmction  with  waves 
bdow  fc  may  produce  ag^cant  CX*'  heating  [59].  Furthermore,  sudi 
iKmresoiuuit  energizatkm  may  wdl  gwe  hi^  enough  energies  to  CTqrlain  some 
km  outflow  fitmi  the  ionoqrhm  [61].  Neverthdess,  Ball  and  Andre  [24,  59] 
conduded  that  cydotrtm  resonance  currently  provides  the  best  explanation  for 
the  miqtnfty  of  observed  km  corac  events  (tl^  are  gmerated  either  locally  or 
atvtr  an  extended  r^km  in  space)  at  altitudes  above  a  few  thtnisand  kilometers 
in  the  aurmal  and  cuq>/ddt  r^kms. 

5.  HEATING  AT  HKHffiR  HARMONICS 

Resonant  heating  of  kms  may  occur  also  at  higher  harmonics  of  the  ion 
g^mbecpuncy  (f==  N/c,  2,3,4...).  As  in  the  case  of  subhamxnuc  heating. 
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energizatkm  at  higgler  hannonics  requires  finite  wavdength  etoctric  field 
fluctuations.  Heating  at  hi^ier  hannonics  is  important,  e  g.,  in  oiergization  of 
Tduunak  fiisitm  plasmas  with  radio  flequency  waves.  Here  waves  are  aitifidally 
launched  on  the  ouigiMtostniic  branch  fiom  the  low  magnetic  field  side  of  the 
plasma.  In  the  nugnetosphere,  wave  al^mption  at  the  second  hamKmic  may  be 
imptxiant,  e.g.,  i^ien  waves  are  propagating  down  a  magnetic  fidd  line,  and  thus 
into  a  strtmga’  magnetic  fidd  [62].  In  this  situation,  the  fi«quaicy  of  a 
downgoing  wave  will  first  mat(^  the  second  harmonic  before  readii^  the 
fimdamental  gyrofietpiency.  Indeed,  in  some  cases  the  downgoing  wave  may  be 
reflected,  e.g.,  at  the  ion-i<ni  hybrid  firequency,  before  reaclnng  the 
gyroflequerwy. 

Observ^  structures  at  higher  hannonics  of  the  proton  gyrofi-equeocy  (N  = 
S  -  10)  might  be  interpreted  as  evidmue  of  absorption  of  broadband  waves.  Ion 
harmonic  structures  (lower  q)ectral  densities  near  multiples  of  fc)  detected  by  the 
S3-3  satdlite  have  been  suggested  as  evidoice  of  energy  transfer  to  the 
smBihaneously  observed  itm  conics  [63],  However,  ion  conics  may  generate 
emissions  between  lugher  hannonics  of  the  ion  ^rofiequency,  induding  waves 
around  the  lower  hybrid  fi'equency  [64].  Here  the  nuqor  source  of  firee  energy  is 
the  positive  slqpe  perpendicular  to  the  geomagnetic  fidd  in  the  itm  vdo^ 
distributitm  Hals,  harmonic  structures  may  also  indicate  wave  gmoation  rather 
than  wave  danqnng.  Structures  at  uppo-  harmonics  of  the  proton  gyrofirequency, 
above  the  lower  favtnkl  firequency,  have  also  been  obsoved  by  the  sounding 
rodtet  MARIE  [86].  Again  these  structures  are  not  necessarily  due  to  wave 
danqmig.  The  peaks  in  the  qrectrum  can  be  H*^  Bernstein  modes,  generated  by, 
e.g.,  parametric  decay  of  lower  l^brid  waves. 

For  suitaUe  conditions,  wave  absorption  at  higher  harmonics  might  be 
significant.  However,  fix  the  broaband  ^>ectra,  which  very  <rflen  are  associated 
with  ion  conics  (see  e.g..  Figure  1),  the  ^ve  intensity  is  much  hi^io’  at  the 
fundamental  km  gyrofiequency  than  at  higher  harmonics.  Thus,  at  high  altitudes, 
ion  cydotron  resonance  herding  at  the  gyrofiequency  usually  is  the  more 
inqxxtant  medumism. 

6.  HEATING  AT  LOW  ALTITUDES  BY  HISS  AND  LOWER  HYBRID 
TURBULENCE 

At  altitudes  bdow  a  few  thousand  kilometers,  cyclotron  resonance  heating 
by  broadband  waves  with  peaks  at  voy  low  fiequendes  m^  be  less  inqmitant. 
At  these  altitudes  the  q>ecttal  density  at  the  Ic^  ion  gyrofirequency  may  be 
rather  low.  This  is  dn^  due  to  the  rapid  decrease  of  spectral  density  with 
firequency  of  the  common  broadband  q)ectra  (see  e.g.  Figure  1).  At  low 
altitiides,  fliere  m«y  still  be  significant  ion  heating  at  the  fundamoital  ion 
gyrofiequen^  by  Inoadband  waves  \^ch  peak  at  rather  high  fiequendes. 
Observdkms  by  the  AUREOL  3  satdlite  at  altitudes  betweoi  one  and  two 
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thousand  Idltmieters  indicate  interaction  between  ions  and  broadband  hiss  waves 

[66] .  These  waves  sometimes  show  a  sharp  lower  cutoff  near  the  proton 
gyrofrequency  at  a  few  hundred  Hz.  Details  of  the  emissions  near  this  cutoff  has 
b^  suggested  as  evidence  of  cooveraon  of  wave  energy  from  one  wave  mode 
to  another  before  ion  heating  occurs,  as  discussed  in  reports  by  Le  Qudau  et  al. 

[67]  and  Le  Qudau  and  Roux  [68],  It  was  noted  by  Johnson  et  al.  [69-71]  that 
the  concentration  of  the  minority  ion  spedes  must  be  low  (a  few  percent)  for  the 
nwdianism  discussed  in  these  rq>orts  to  be  significant.  Furthemore,  at  these 
altitudes  the  gradioit  scale  lengths  may  be  of  the  order  of  the  wavdength,  and 
thus,  e.g,.  boundary  conditions  can  berome  important  in  estimates  of  the  wave 
abs(»ption  rate. 

More  narrow-banded  waves  might,  at  least  in  principle,  also  energize  ions. 
There  are  some  rocket  observations  of  O'*'  cyclotron  harmonic  waves  together 
with  perpendiculaiiy  heated  ions  [3].  However,  the  waves  do  not  seem  to  be 
intense  enough  to  e3q)lain  foe  observed  ion  heating. 

It  seems  likely  that  emissions  around  and  above  the  Iowa-  hybrid  frequent^ 
are  inqx>rtant  for  ion  energization  at  low  altitudes,  particularly  in  the  region  of 
discrete  aurmas.  Several  years  ago^  Chang  and  Coppi  [72]  su^ested  that  lower 
hybrid  turbulence  could  be  the  prime  candidate  for  tte  transverse  accderation  of 
ions  in  the  high-latitude  topside  iono^here,  particularly  in  the  r^on  of  discrete 
auixvas.  Since  precipitating  energetic  dectron  beuns  are  assodated  with 
discrete  auroral  diqdays,  it  was  argued  that  these  dectron  beams  could  provide 
the  free  eneitgy  fi»  the  production  of  intense  Iowa*  hybrid  waves  that  are 
observed  there.  The  ctHnmonly  observed  geometry  of  disoete  auroras  is  that  the 
arcs  are  extended  in  the  looghudinal  direction  but  narrowly  confined  in  the 
latitudinal  direction.  This  indicates  that  foe  r^on  containing  the  energetic 
precipitating  dectrcms  above  a  discrete  auroral  arc  (the  suprauroral  r^on)  must 
also  be  similariy  confined  in  the  latitudinal  direction.  The  dmsity  whW  foe 
suprauroral  r^on  is  generally  less  than  tiiat  of  its  immediate  nei^iboihoods  with 
af^nodable  perpendicular  (outward)  and  parallel  (downward)  dentity  gradients. 
Tlie  ixropagation  character^tics  of  the  rdevant  plasma  waves  that  are  excited  by 
the  energetic  electron  beam  at  high  altitudes  depend  oudally  on  the 
inhmnogeneous  density  structure  of  the  density  depleted  suprauroral  region.  It 
has  been  shown  by  Majggs  [73],  Maggs  and  Lt^o  [74],  and  Retterer  et  al.  [75], 
using  ray-traditg  calculafoms,  that  for  typcal  plasma  parameters  in  foe 
suprauroral  r^ion,  convective  modes  at  fi’equendes  much  higher  than  the  lower 
hybrid  fiequency  will  generally  propagate  out  of  the  r^on  in  the  latitudinal 
direction  shortly  after  they  are  graereted.  Thus,  only  those  modes  dose  to  the 
lower  hybrid  fiequency,  vdnch  have  group  vdodties  neariy  paralld  to  the 
magnetic  fidd,  can  continue  to  propagate  and  grow  along  the  fidd-aligned 
direction.  These  lower  hylmd  mo^  can  stay  whlun  foe  suprauroral  region  and 
grow  to  very  large  artqrlifodes  leading  to  the  nonlinear  coupling  and  cascading  of 
modes  into  a  Inoadband  spectrum  qranning  frcmi  the  origirud  lineariy  unstable 
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wsvdengths  to  shorter  and  shorter  wavelengths  at  lower  and  lower  altitudes 
[76,77]  Figure  11. 

It  was  shown  by  Retterer,  Oiang  and  Ja^>a8e  (RCJ)  [76]  that  such  type 
of  ncHilinear  coupling  process  may  be  understood  in  terms  of  the  phoiomenon  of 
nKxhilatkmal  instability  (due  to  the  sdf-consistently  genoated  low  frequency 
modes)  based  on  the  warm  fluid  approximation  for  the  ions  and  dectrons.  In 
addition  to  the  nonlinear  coiq)lii%  process  of  wave-wave  interactions, 
ccmomnitant  wave-particle  intoactimis  control  the  evolution  of  the  Iowa-  l^rbrid 
turbulence  in  the  lower  altitude  magnetoqrfxre  and  topside  iono^tme.  RCJ 
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Figure  11.  Typical  wave  qwctia  from  a  one-diineasional  lower  hybrid 
PIC  shnnlatkn  at  few  diflfennt  times  with  t4  >  t3  >t2  >ti.  Itisseenthat 
the  wave  energy  cascades  cmitinually  from  the  k-vahies  of  the 
linearly  unstable  modes  toward  lar^  k-vahies  indicating  the 
casca^g  process  toward  dxnter  wavdength  lower  hybrid  waves. 
From  Rettm  et  al.,  1S>86,  [76]. 
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suggested  that  the  dominant  wave-particle  interaction  process  could  be  modeled 
by  the  consideration  of  instantaneous  Landau  damping  of  the  waves  by  the 
various  particle  species  (i.e.,  the  beam  electrons,  the  background  ambient  ions 
and  dectrons)  which  at  the  same  time  are  scattered  by  the  broadband  lower 
hybrid  waves  through  resonant  interactions. 

Figure  12  depicts  the  evolution  of  the  amplitude  of  the  dectric  field  in  real 
space  and  time  based  on  the  RCJ  model  [78,  79].  At  the  beginning,  the 
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Figure  12.  Evohitiaa  of  lower  hylmd  caviton  turbulence  in  real  space  and 
time  based  on  die  RCJ  [76]  dieory.  Shading  indicates  die  strengdi  of  the 
electric  field  amplitude.  The  hcmzontal  scale  is  approximately  100  D^ye 
lengdis  and  die  vettical  scade  is  approximately  100  ion  plasma  periods. 
Nodoe  how  die  initial  nearly  unifimn  wave  anplitude  intensity  ewdves  into 
intense  individuad  localized  solitary  structures  or  broadband  spiky  (cavitcn) 
turbulence.  From  Chang,  1993,  [79]  and  Retterer  et  ad.,  1993,  [78]. 

anplhude  of  the  electric  field  grows  unifonnly.  Eventually,  when  wave  cascade 

sets  in,  the  electric  field  intensity  begins  to  concentrate  ^periodically  and  in 
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narrow  qMtial  r^<Mis.  These  localized  structures  (cavhons,  because  they 
conespoiid  to  localized  density  depledcnts  [76,  80])  are  seen  to  propagate  at 
q)eeds  much  less  than  the  ambient  ion  thennal  velodty  and  should  appeu  as 
staticHiary  structures  if  detected  by  in  situ  electric  field  measurements. 
Smnetimes  these  structures  are  densely  mixed  forming  a  ^nky  turbulence  with 
Ixoat&and  waveki^ths  modes.  At  low  altitudes  (topside  ionoq)here;  MARIE 
and  TOPAZ  attitudes),  file  individual  localized  structures  that  are  dominated  by 
slxvt  wavder^  lower  hybrid  modes  can  initiate  the  transverse  energizatiou  of 
ionoqrfieric  ions.  At  mid-altitudes  in  the  suprauroral  r^on  (e.g.,  ISIS-2  or 
Frqa  altitudes),  the  itms  energized  by  the  cavhons  or  other  preheating  processes 
can  be  heated  confimially  1^  the  broatfiiand  (spiky)  lower  hybrid  turbulence  to 
hitler  energies.  These  broadband  lower  hybrid  waves  can  also  energize  the 
ambient  electrons  in  (both  the  upward  and  downward)  fidd  aligned  directions. 
The  individual  localized  structures  at  mid-altitudes,  on  the  other  hand,  may  not 
fiid  suflSdentfy  cool  ions  to  resonate  whh.  They  can,  neverthdess,  interact  with 
the  ambient  dectrons  if  the  resorumt  contorts  are  met,  and  produce 
counterstreaming  bursts  of  energized  dectrons. 

Recent  rodeet  experiments  [4,  6,  7]  have  obtained  high-time  resolution 
wave^Mutide  data  in  tte  topside  iono^hm-e,  that  are  in  remarkable  agreement 
with  tte  predictirms  of  the  Tories  of  Chang  and  Coppi  [1981]  and  RCJ  [1986]. 
Details  of  the  results  of  these  innovattve  expoiments  on  the  observation  of  Iowa- 
hybrid  cavitems  and  assodated  transvm^  ion  accdoatimi  are  given  in  separate 
reviews  [81, 82]. 

In  addition,  the  difiuskm  themy  of  broadband  q;)iky  turbulence  of  an 
ensemMe  of  cavhmis  (or,  shi^,  broa^iand  Iowa-  faylnid  turbulence)  was  fiiund 
to  be  adequate  to  esqilabn  the  observed  ion  crmics  the  S3-3  sat^e  and  the 
MARIE  rodeet  [75].  A  difiiision  theory  based  on  the  idea  of  comUned 
preheating  by  c^lotron  resonance  and  continued  oiergizatirm  by  broarfliand 
lower  hybrid  turtnilence  has  also  given  a  plausible  eiqilanation  of  the  recently 
observed  transversdy  energized  ionosfrfieric  ion  (0+)  events  by  the  Frcja 
satdlite.  [See  discussitms  given  in  Section  10.] 

Findly,  it  is  to  be  noted  that  lowor  hybrid  waves  can  be  genoated  by  other 
mechanisms  sudi  as  ion  ring  distributions  [127],  ion  drifts  [128],  and  vdodty 
shear  [129].  In  fiict,  ootain  events  of  transverse  ion  acedoation  obsoved  by 
S3-3,  DE  and  Fr^  satdlites  have  been  identified  with  the  sunultaneous 
occurrence  pre^ntating  protm  and  hdium  rings  [127, 130]. 

7.  SOME  ADDITIONAL  EXAMPLES  OF  ION  HEATING  MECHANISMS 

7. 1  Electrostatic  Ion  Cydotron  Waves 

Sevml  possible  itm  heating  mechanisms  in  addition  to  those  discussed 
above  have  been  considered  by  various  authors.  One  mechanism  that  has 


received  considendde  atteodcHi  is  heating  by  electrostatic  ion  waves  above 
muh^les  of  the  ion  gyrofiequaicy.  These  waves  noay  be  generated  by 
simultmeously  observed  ion  beams,  or  possibly  1^  drifting  dectrons 

[83-87].  Both  analytical  calculatimis  [e.g.,  88]  and  siinulations  [89;  and 
references  therein],  show  that  electrostatic  km  cydotrtm  waves  may  cause 
perpendicular  ion  heating.  However,  there  is  little  observational  evidence  that 
this  energizatkm  process  is  of  any  miyor  importance  in  the  terrestrial 
magneto^here.  At  high  altitudes,  th^  waves  are  usually  observed  together 
with  ion  beams,  rather  than  ion  conics. 

There  are  rodcet  observations  of  oxygen  cyclotron  harmonic  waves  and 
perpendiculaily  heated  O'*’  kms  [3],  but  th^  waves  do  not  seem  to  be  intense 
enough  to  cause  die  observed  ion  heating.  However,  km  beams  may  have  large 
perpendicular  tenqieratures  (as  comqiared  with  mean  energies  in  the  Iowa* 
kmo^here),  larger  transverse  than  paralld  temperature  [e.g.,  90]  and  also  a 
vague  conic  duyie  [e.g.,  86].  This  noay  indicate  some  perpendicular  ion  heating 
by  broadband  low  firequency  waves.  Furthermore,  in  many  situations  the  most 
plausiUe  candidates  may  be  the  acoustic  waves  as  investigated  by  Bergmann  et 
al.  [91].  This  is  further  discussed  bdow,  in  connection  with  interaction  between 
i<m  beams. 

7.2  Trapfrag  by  Intense  Cohoent  Waw»  and  Stochastic  Accderation 

When  one  wave  mode  donunates  the  wave  fidd,  ion  tnqiping  in  this  mode 
may  occur.  In  general,  the  ion  trapping  mints  will  limit  the  energy  achievable  by 
the  kms.  However,  as  discussed  in  a  review  by  Lysak  [22],  the  interactions  of 
the  kms  with  higher  harmoiucs  or  with  other  wave  modes  may  allow  ions  to 
junqi  frmn  one  trapi^  oiiA  to  another.  This  can  lead  to  an  irreversilde 
increase  in  the  oiergy  of  the  ions.  For  this  modd  to  be  meaningfiil,  the  wave 
fidd  must  be  coher^  and  sufiBcmntly  strtmg  [92-94].  The  observations  of 
strong,  narrowbanded  dectrostatic  ion  cydotron  waves  within  the  auroral 
accderation  r^on  by  Mozer  et  al.  [95]  motivated  Lysak  et  al.  [96], 
Papadopoulos  et  al.,  [97],  and  Singh  et  al.,  [98,  99]  to  consider  stochastic  ion 
heating.  However,  it  is  not  dear  that  coherent  and  intense  waves  are  sufiSdently 
cmnmon  in  the  terrestrial  magneto^here  for  this  process  to  be  of  any  miyor 
inymitaime. 

7.3  Intoaction  Bdweeo  Ion  Beams 

For  several  of  the  ion  heating  mechanians  discussed  in  this  report,  we  have 
not  considered  the  wave  generation.  In  the  auroral  zone,  the  situation  is  rather 
conyilex,  but  <me  possible  generation  mechanism  is  interaction  between  ion 
beam.  A  static  potential  drop  accderatii%  ions  upward  will  give  all  ions  the 
same  oiergy,  and  thus  the  li^er  ions  will  get  higher  vdodties  [100;  and 
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reftreoces  therein].  Ion  beams  obsenred  at  higher  ahhudes  have  much  larger 
tcnqMralures  than  the  soince  ions  [e.g.,  90].  Thus,  these  ions  are  heated  as  wdl 
as  aooderaled  in  the  potential  structure.  One  posaMe  source  of  wave  generation 
and  sidMoquent  transverse  km  energiation  are  instaUlhies  caused  by  the  rdative 
drift  betw^  these  ion  beams  [87,  101*104].  Simulations  show  that  the 
lesidting  wave  enusskms  can  cause  scmie  eoei^zation  both  in  die  paralld  and  the 
petpeadicuiar  directkm  [105-108;  and  reftrences  therein].  Scmie  of  the 
perpendicular  heating  may  be  relevaitt  ft>r  our  understimding  of  the  i<m 
distributions  usually  called  ccmics.  Sevoal  other  wave  generation  mechanisms 
are  posskde  on  dieaefidd  lines.  The  pMsibilities  indude  interaction  between  km 
beam  and  a  statkmaiy  ion  badcground,  eiectnm  beams,  and  varkms  gradients. 
Thus,  again,  observaticms  are  needed  to  identify  the  important  wave  generation 
medumiam(s). 

7.4  Heating  of  Hdhun  in  the  Equattxial  Kane 

In  the  equattnial  |dane  near  geostationary  mbit  both  wave  endsaiims  and 
their  energy  source,  as  wdl  as  heated  icms,  m^  occur  in  the  same  r^cm. 
Anisotropic,  hot  (MXiton  distrSmticms  [109,  1 10]  mi^  generate  dectnmaagnedc 
km  cydotron  waves  bdow  the  inotrm  gyroftequi^  [109,  110-115].  As 
disaiMed  bdow,  such  waves  may  pr(^)agate  down  to  the  central  fdasma  sheet. 
Ifowever,  the  waves  mi^  also  rather  locally  heat  He'*'  ions  [1 12].  We  do  not 
egqiect  ftiese  hdium  distributions  to  have  a  crmical  shape,  but  we  indude  them  as 
an  examine  of  transverse  ion  energization. 

In  the  ftdlowmg  section  we  discuss  how  waves  may  propagate  fron  the 
equatorial  plane  and  cause  ion  heating  at  lower  altitudes. 

7.5  Wave  Generation  and  Propagatkm  in  the  Central  Plasma  Sheet 

BroadMnd  waves  and  km  conics  are  often  observed  in  the  r^on  of 
discrete  auroras  and  in  the  cu^deft  r^km.  Here  several  local  energy  sources 
sudi  as  sharp  gradients  and  drifting  partides  can  possiUy  genmate  some  of  the 
waves.  However,  broadband  waves  associated  with  ion  conics  are  also  commem 
in  the  central  plasma  sheet  [41].  As  an  emmqrfe  of  wave  generatkm  and 
transport  associated  with  km  conks  we  discuss  srnne  recent  investigations  of  the 
central  plasma  sheet.  In  this  r^<m,  equatenward  of  the  r^km  of  disoete 
auroras,  there  are  no  obvious  local  energy  sources  that  can  power  the  broadband 
waves.  This  led  Johnstm  et  al.  [69]  to  suggest  that  the  waves  are  generated  by 
niaataofic  km  distributions  in  the  ecpiatorial  fdane.  It  is  wdl  known  ftiat  left- 
Inmd-poiaiized  waves  bdow  the  protem  gyrofiequency  can  be  generated  in  the 
equatenial  ibme.  Hoe  a  mixture  of  cool  and  n^vdy  dense  {dasma  with  hot 
and  anisotropic  |dasma  is  rather  common.  The  hot  ion  oompoaaat  often  has  a 
hqkcr  temperatine  perpendicular  than  paraUd  to  the  amUent  magnetic  fidd 


[109],  <Mr  a  significaitt  loss  cone  [110].  Sudi  a  mixture  of  cool  itms  and  a  hot, 
anisotropic  km  population  is  fin^wrable  for  the  amfdification  of  km  cyclotrcm 
waves  below  the  prot<m  gyrofrecpiency,  as  discussed  by,  for  aaaapk,  Cornwall 
[111].  These  waves  may  then  {xopagate  down  the  geomagnetic  fidd  lines.  As 
discussed  by  Johnson  et  al.  [69],  conversion  of  wave  energy  from  one  wave 
nmde  to  anodier  m^  be  importmit  as  the  ^ves  pr<^>agste  downward. 

Similar  km  cydotnm  waves  in  the  plasnu^Muise  r^km  have  been  studied 
by  Home  and  Thome  [62],  who  used  ray  tracing  to  study  the  wave  prc^Mgation. 
lliey  emphasized  km  cyclotron  danq>^  at  the  second  harmonic  of  the  O'*' 
gyrdfrequen^.  The  pafo-intcgrated  d>sorpti<m  they  conqmted  was  used  to 
estimate  the  km  heating  qualitativdy,  but  th^  method  did  not  allow  quantitative 
conqNuisons  based  on  ohierved  spectral  denaties. 

In  a  recent  study,  Rbmmiaric  and  Andr6  [1 16]  discussed  the  convection  of 
km  cydotron  waves  fram  the  equatorial  i^ane  to  lower  altitudes  where  heavy 
kms  are  heated,  hi  this  investigatkm,  a  moddofthe  equatorial  plasma  is  used  to 
estimate  linear  wave  growth  and  to  modd  the  distri^rm  of  wave  energy  in 
phase  qmoe.  Phase  qmoe  methods  devdqied  by  Rimnmaric  and  Larsson  [117] 
are  then  used  to  nug)  qiectral  densities  along  the  magnetic  fidd  lines.  The  ray 
traciQg  mediods  us^  by  Rtenmark  and  Andrd  [116]  are  rather  conventional. 
Ifowever,  it  is  emphadzed  that  the  are  curves  in  fdiase  qiace.  The 
conservation  of  the  so  called  wave  distrfoutkm  fonction  akx^  rays  is  then  used 
to  andyze  the  transput  of  wave  energy.  This  method  makes  it  possible  to  rdate 
the  spectral  density  observed  in  die  iofr^ieatmg  rq(km  to  measured  spectra  in  the 
equatorial  |dne.  The  conduskm  obtained  by  Rdonmaric  and  Andrd  [1 16]  is  that 
wave  energy  oonvecting  down  the  fidd  lines  in  the  central  plasma  sheet  can 
significant^  contrfoute  to  the  waves  observed  in  r^kms  of  ion  energization. 

More  recently,  Jolmson  et  al.  [70,  71]  demonstrated  that  magnetosonic 
waves  could  also  be  ducted  fixnn  the  equatorial  r^on  to  auroral  hdghts  by 
peqpendiailar  density  gradients. 

We  end  the  discusskm  of  wave  generatirm  and  propagation  here.  In  the 
next  section  we  briefly  oondder  paialld  km  accderation  that  is  often  rdated  to 
perpeuficular  km  eneirgizatkm 

8.  PARALLEL  ACCELERATION  ASSOCIATED  WITH  ION  HEATING 

Ions  ftxming  so  called  km  ctmics  often  display  diaracteristics  of  not  <mly 
perpendicular  but  also  of  paraUd  ener^zation.  This  signature  consists  of  fidd 
afigament  of  die  kms  at  energies  just  above  a  lower  energy  cutoff  bdow  winch 
essential^  i»  kms  ate  observed  (e.g,  2).  At  itxxeasing  energy,  the  fidd> 

alined  distribulion  gradually  beomies  a  ccmi^  one  whose  qmx  ai^  widens 
with  increasing  energy.  These  distrfoutions  are  sometimes  called  bimodal,  or 
bowi-diaped,  or  lifted  km  conics.  The  first  discussion  df  these  km  distributions 
suggested  a  two  stage  process  [1 18].  Ions  may  first  be  transversdy  energized. 
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nd  then  acedonted  throng  a  I>araUd  potential  drop.  It  has  also  been 
suggested  thtt  a  panifld  electric  field  that  may  be  deconqiosed  into  a  fluctuatiog 
partandastaticpartmaygivesiimlar  effects  [15, 119].  When£j|  is  changing  <» 
a  suittMe  timescale,  ions  nu^  experknoe  essentially  an  average  downward 
electric  fidd.  The  dectrems,  however,  have  mudi  larger  vdodties.  These 
particles  may  sti^  in  the  re^on  of  interest  during  less  than  one  wave  cycle,  and 
may  be  accdeiided  either  upward  or  downward  (given  a  large  enougli 
fbetuming  part  of  £||).  Such  thne>varying  £||  is  (me  possfide  exidanatkm  of 
observations  (rf'krns  and  electrons  moving  iqnvardtoget^.  However,  we  note 
that  the  km  aoederatkm  by  a  panOd  dectric  fidd  does  not  give  perpendicular 
encfgizatkm.  Thus,  the  cemied  diqm  (distorted  by  paraUd  accderati<m)  has  to 
be  expidned  sqmratdy. 

It  is  not  necessary  to  assume  paraOd  acederatitm  of  the  whole  distributiem 
to  cause  lifted  conics.  Rather,  the  low  energy  part  of  the  conic  may  be  removed 
on  tile  way  to  the  point  of  observatkm.  This  nu^  ()ccur  via  vdkmity  filtering 
[120].  In  this  scenario,  the  km  dtstrSnitkm  is  transversely  bulk  heated  in  some 
finiter^km.  TUs  process  is  ftdowed  by  essentially  adiabatic  convective  flow  to 
the  observation  kxxttkm.  Here,  the  latter  stage  may  contribute  a  velocity  filtering 
effect.  Again,  the  perpendicular  heatiqg  and  the  liftiaig  of  the  distrifeutkm  in 
paralld  vdocity  are  essentially  separate  processes. 

Some  km  heatup  medhanisms  may  be  deserfted  as  a  randcmi  walk  in 
perpendicular  vdocity.  Su(di  a  process  together  with  iqpward  motion  in  a 
magnetic  fidd  with  decreasing  nu^pntudei,  can  also  give  itm  conics  that  are  lifted 
in  pardM  vdocity  [121].  Note  tiiat  <mty  perpendicular  energization  acting  rm 
the  idKde  tbtritxaion  is  needed  in  this  case.  At  first  it  may  seem  strange  that  a 
process  that  accelerates  kms  only  perpendiculariy  to  the  magmeifi  fidd  can 
prckhice  a  (hstmet  low  vdocity  cut^  in  peraDd  km  vdodty.  However,  as 
diacusaed  by  Temerin  [121],  heating  »  a  stodiastic  process.  That  is,  tiiou^  the 
healed  km  distributkm  always  gain  energy,  individud  kms  nuty  either  gain  or  lose 
energy  at  each  tme  atq>  in  the  random  walk.  Ions  with  wwall  perpendicular 
vdocities  observed  by  a  satellite  may  have  reedied  tins  qmcecraft  ^  moving  up 
the  fidd  fine  through  a  long  heating  r^km.  During  the  time  tlm  kms  were 
moving  up  the  field  fine,  tii^  may  have  gained  and  kmt  perpendkular  energy 
several  times.  During  the  times  they  had  sigmficant  perpendicular  energy,  they 
gained  paraOd  vdocky  due  to  the  magnetic  mirror  fi:^.  Thus,  it  is  extremely 
improbiUe  that  an  km  dxwld  never  have  gained  aity  s^nificant  perpendicular 
energy  ifthere  are  a  large  mutiber  of  inleracticxn  in  the  random  wafic.  Hence,  all 
kms  dwidd  have  gained  some  parafld  velocity.  As  point  out  by  Chang  etal.  [41] 
and  Retlerer  et  d.  [43]  and  discussed  earlier,  cydotron  rescmance  hft«*w»g  by 
broadbmd  waves  is  one  (^  Ae  processes  that  may  be  described  as  a  randixn  walk 
in  perpcndicidar  vdocity.  Furtiiernxxe,  heatog  over  an  extended  altitude  range 
is  naa^  fly  Ititing  m  parafld  velocity  to  bectmae  obvious.  Thus,  the  cydotron 
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resonance  hem^  mechanism  is  consistent  with  the  observaixms  of  lifted  km 
conics  at  altitudes  above  a  ftw  thousand  kikmwters  (see  e.g..  Figure  2). 

Recently,  Li  and  Temerin  [122]  si^gested  that  the  ptmdennocive  fi>rce 
derived  from  the  very  low  frequency  portkm  of  the  decthc  fidd  wave  q>ectra 
am  provide  fiekMignedenergtcation  of  kms.  This  idea  inilizes  the  gradi^  of 
the  E^B  drift  in  die  Earth's  diverging  magnetic  field.  The  mechanism  requires 
rather  large  low  frequency  (mudi  bdow  fg)  fluchutting  perpendicular  dectric 
fidds.  A  similar  process  ft>r  parsfld  km  energizatkm  in  tte  Earth's  nonunifbrm 
static  dectrk  and  magnetic  fidds  was  descrftmd  by  Cladis  [1986].  Here  no 
explicit  tune  dcpeodenoe  of  die  fidds  are  introduced  in  the  modd.  Rather,  the 
pveOd  km  aocderatkm  is  propmtkmal  to  the  E^B  drift  and  to  the  rate  of 
turning  of  the  magnetic  fidd  directxm  akmg  the  guiding  center  path  of  the  i<m. 

Another  simiiar  process  was  dmcussed  by  Mauk  [124].  In  this  case  a 
sudden  displacemeot  of  curved  fidd  lines  resulting  from  a  shmt-lived  api^catkm 
of  an  intense  dectric  fidd  may  violate  the  seomd  adiabatic  invariant  of  the  lower 
energy  ions.  Associated  with  das  vidadtm  is  the  condhkm  that  the  particles  that 
gain  the  largest  muldplication  in  their  energies  are  those  at  small  i^di  an^es. 
Depending  on  parameters,  this  mi^  9ve  stgm&ant  paralld  ion  energizatkm. 

In  some  situations  diese  pardld  heating  processes  may  need  to  be 
considered  tt^edier  with  perpendiailar  ion  heating  mechanisms  u^hen  describing 
the  fixmatkm  of  km  conics. 

Some  recent  radar  observatkms  at  altitudes  up  to  1500  km  indicate 
significant  km  outflows,  often  associated  with  auroral  arcs  [125;  and  references 
therein].  These  ions  have  temperatures  of  a  few  thousand  KdKin,  and  field> 
digned  velocities  iq>  to  dmut  one  kilometer  per  second.  The  km  distributions 
have  no  obvious  temperature  anisotropy,  and  thus  are  not  ion  amics.  However, 
these  km  outflows  can  8tq)|dy  particles  to  r^ons  oi  perpendicular  ion  heating  at 
higher  altitudes. 

9.  RECENT  OBSERVATIONS  BY  THE  FREJA  SATELLITE 

The  Fr^  satdlite  was  launched  in  October  1992  and  has  an  apogee  of 
1760kminthenotthenihenB4diere[131].  The  Frga  orbit  with  an  incUnadon  (rf* 
63  is  often  aptnoximately  tangential  to  die  auroral  oval,  making  long  continuous 
observatkms  in  ion  hearing  regions  possible.  Intoise  ion  energizaritm  to 
diaracterisric  energies  ofup  to  about  50-100  eV  (approximatdy  the  same  for  H^ 
and  O^  has  been  obsoved  <m  several  occasiems.  These  ion  conics  are  often 
associated  widi  intense  lower  hybrid  waves  (with  anqilitudes  up  to  and  above 
100  mV^i)  but  the  correbriem  nu^  be  &r  from  perfect.  Usually,  onisskms  at 
low  frequencies  are  also  associated  wfth  the  conics.  The  km  pit^  aigfes  dx>w 
that  most  dftlie  hearing  has  occurred  within  a  few  hundred  km  in  altitude.  One 
Fr^  hearing  event  in  die  eveningside  aurmal  ztme  has  been  studied  in  detail  by 
Andrd  et  al.  [132].  During  this  event,  most  of  the  intense  ion  heating  is  probably 
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cauMd  liy  the  obfoved  tower  hybrid  waves.  These  waves  may  be  caused  by  the 
shmhaneoua^  observed  downgoii^  auroral  electrons  with  energies  of  several 
ImV.  However,  stMoe  additional  energizatkHi  ("preheating")  by  e.g.,  cydotron 
resonance  heatiag  caused  by  waves  observed  around  the  km  gyrofrequendes 
(audh  as  diat  n^gested  by  Chang  et  al.  [133])  seeoos  to  be  needed.  It  should  be 
noted  diat,  in  contrast  to  recent  observatkms  by  soundii^  rockets,  many  Frqa 
km  heating  events  are  not  associated  with  indivkhud  cavhons  [134].  This  is 
consistent  with  the  discusskms  given  in  Section  4.  The  process  ^  modulatkmal 
mstability  can  produce  broadband  wavdength  tower  hybrid  waves  with  or 
without  distinct  indhndual  localized  structures  (or  cavhons).  TUs  is  mie 
irrespective  of  the  generatkm  mechanisms.  As  asserted  earlier,  totmdband  tower 
hybrid  waves  and  transverse  km  aocdentfkm  has  also  been  observed  in 
coiymctkm  with  predpitating  km  rings  by  Frqa[  130]. 

10.  DISCUSSION  AND  CONCLUSIONS 

Ion  conics  are  often  detected,  but  sumihaneous  observatitms  of  waves  and 
local  transverse  heating  Oon  diatributkms  which  peak  at  90*^  pitch  an^)  by 
sateBites  at  Ugh  altitudes  are  not  common  [126].  Such  distributions  are  only 
occaatonally  observed  at  altitudes  above  a  few  thousand  kilometers  [13], 
ahhoi^  several  years  of  satellite  data  are  avaihtole.  Some  exanqdes  of 
esaentiafly  90**  conics  have  been  ftwnd  by  rodcets  at  tower  altitudes  [4,  5,  6, 1 1, 
12;  82].  However,  mediaiBsms  such  as  electrwn^netic  km  cydotron  resonance 
operate  over  thwisands  of  Idtometers  in  altitude,  gradually  heating  the  ions, 
liais  it  is  not  surprising  that  only  a  few  locally  heated  ion  disritoutions  have  been 
ftamd  at  U^  altitudes. 

It  seems  likely  that  at  least  a  few  diffiaect  ion  heating  medumisms  (wave 
modes)  are  important  in  die  terrestrial  mi^netomdiere.  Future  tests  of  these 
mecharisms  dmild  involve  events  vriiere  mmre  tl^  ooe  heated  itm  q>edes  is 
observed  at  a  specific  time.  The  hearing  correspondii^  to  a  certain  mechanism 
and  ton  mass  can  be  estimated  from  an  observed  wave  plectrum.  However,  the 
calculation  often  involves  several  assumptions  (e.g.,  about  wavdei^  and  wave 
polarization).  Events  with  more  than  one  km  species  should  provide  even  betto* 
tests  of  the  dwoiies.  Statistical  studies  (rf*  km  conic  events  i^ere  both  wave  and 
partide  data  are  avaihMe  would  alao  help  to  clarify  which  mechanisms  are 
important  in  practice.  AH  these  observarional  studies  require  h^  resdution 
data  fioni,e.g.,  die  auroral  zone.  One  source  of  such  data  is  die  Fr^  satellite, 
widimiiqiogee  of  1760  km.  Thus,  tUs  ipaoecraft  can  provide  a  large  amount  of 
^ta  fiom  a  region  where  H^,  He'*',  and  ions  ve  common.  Odier  interestmg 
data  w9  be  provided  by  e.g.,  die  FAST  satellite,  planned  to  be  launched  durug 
1$194,  and  by  soundmg  resets.  All  dieae  data  shoidd  be  used  not  only  to 
investigate  detads  of  the  km  heatmg  mechanisms,  but  also  to  increne  our 
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understanding  of  large  scale  phenomena  such  as  ionoqirfiere-inagnetospha’e 
coupling  and  nonthermal  esctqpe  of  planetary  atmo^ho-es. 

We  have  discussed  a  number  of  i<m  heating  mechanisms  applicable  in  the 
terrestrial  magnetoq)here.  At  low  altitudes,  energization  by  turbuloice  around 
the  lower  hybrid  frequency  seem  to  be  one  iix^HMtant  heating  mechanism.  Some 
kms  that  have  been  perpoidiculaily  heated  at  low  altitudes,  may  be  observed  at 
higher  altitudes  as  trowing  "clas^cal  comes”  with  their  vertices  located  at  the 
migin  in  vdodty  space.  At  altitudes  above  a  few  thousand  kilmneters,  cyclotron 
resmumce  heati^  by  Ivoadband  waves  armmd  the  ion  gyrofrequency  seem  to  be 
the  most  inq)ortant  process  for  transverse  ion  ener^zation.  This  mechanism 
provides  bulk  heating  of  the  ion  distributions.  Consequently,  the  conics  resulting 
frmn  heatup  over  an  extended  altitude  r^on  will  generally  be  lifted  in  paralld 
vdocity.  All  these  ion  heating  mechanisms  contribute  to  the  outflow  of 
icmospheric  plasma  into  the  nu^petosphere. 
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